Zeno and anti-Zeno polarization control of spin-ensembles by induced dephasing 
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We experimentally and theoretically demonstrate the purity (polarization) control of qubits entangled 
with multiple spins, using induced dephasing in nuclear magnetic resonance (NMR) setups to simulate 
repeated quantum measurements. We show that one may steer the qubit ensemble towards a quasi- 
equilibrium state of certain purity, by choosing suitable time intervals between dephasing operations. 
These results demonstrate that repeated dephasing at intervals associated with the anti-Zeno regime lead 
to ensemble purification, whereas those associated with the Zeno regime lead to ensemble mixing. 
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Introduction. — The ability to understand and manipu- 
late the dynamics of "open" quantum systems, i.e. sys- 
tems that interact with their environment ("bath"), is a ma- 
jor challenge for fundamental quantum physics and a pre- 
requisite for novel applications such as quantum heat en- 
gines lHll, quantum information storage and retrieval 10] 
and precision measurements fs]]. This is particularly true 
of information-carrying spin-1 /2 particles, known as quan- 
tum bits (qubits), coupled to spin-1/2 particles of other 
species that constitute the bath. Such systems are usually 
controllable by coherent fields 10-01 ■ Here we address their 
manipulations by incoherent, random fields that dephase 
the system and thereby mimic quantum non-demolition 
(QND) measurements lH Such manipulation is con- 
ceptually intriguing: whereas QND measurements leave a 
closed system intact, they can affect an open system, by 
destroying its correlations (coherences) with the bath. As 
recently predicted for qubits coupled to thermal oscillator 
baths, such measurements can steer the qubit ensemble, to- 
wards either higher or lower purity ("cooling" or "heating") 
ifioll : Namely, the qubit does not retain its state as the mea- 
surements accumulate, but rather converges to an asymp- 
totic steady state. In the Quantum Zeno (QZE) regime, 
highly frequent measurements raise the asymptotic excita- 
tion and entropy (mixing) of the qubit. This reflects the 
hitherto unnoticed fact that QZE dynamics equalize the 
bath-induced upward and downward transition rates, in the 
qubit. By contrast, less frequent measurements conform- 
ing to the anti-Zeno (AZE) regime ficll, predominantly 
enhance downward transitions (relaxation to the ground 
state) and thus purify ("cool") the qubit. These predicted 
measurement-induced changes of the equilibrium go be- 
yond previous studies that focused on transition-rate (re- 
laxation) slowdown in the QZE regime and its speedup in 
the AZE regime |@] and have been experimentally verified 
ill. 



setups Jt]], and field-driven quantum dots 1 12|. Since all 
the / spins have the same energy levels, such spin-baths 
are spectrally degenerate, as oposed to the broad spectrum 
of oscillator baths. The resulting qubit-bath dynamics is 
therefore different for the two scenarios and hence we ask: 
do the equilibrium changes predicted in Ref. ifloll hold for 
both scenarios? 

In this work we demonstrate that they do, despite their 
differences: we experimentally lower or raise the purity of 
the system and bath spins via frequent induced dephasings 
that simulate QND energy measurements |0] respectively, 
by timing the dephasing intervals to be in the QZE (evo- 
lution slowdown) or AZE (evolution speedup) regime |@] 
respectively. Remarkably, repeated dephasings at intervals 
conforming to the AZE lead to highly-effective polariza- 
tion exchange that can overcome even large frequency de- 
tunings (off-resonant mismatch) of the qubit and bath spins 
and induce polarization transfer that is as large as if it were 
the Hartmann-Hahn resonant transfer ifisll . This novel ef- 
fect, demonstrated by our NMR polarization transfer ex- 
periments, is termed here incoherent resonance, as it stems 
from repeated system-bath correlation erasure (dephasing). 

Model and dynamical regimes. — We assume that the 
qubit-bath system is described by an effective Hamiltonian 
H = Ho + Hsi + HM{t). (1) 

having the following terms: (i) The Hq Hamiltonian ac- 
counts for the coherent evolution of the qubit and the bath; 
under a Zeeman-like interaction with respective Larmor 
frequencies ujs and ojj. (ii) The Hsi term describes the 
coupling between the S and I spins, chosen for simplicity 
to be oriented perpendicular to the Zeeman field, 

Hsi = JY.SUI, (2) 



The present study considers a scenario different from 
that of Ref. ifioll : the interaction of a spin-1/2 or qubit 
system S with N identical spin-1/2 particles / that con- 
stitute its "bath". Such a situation is encountered in NMR 



S'^ and being the a^. Pauli operators for the respective 
species, (iii) The time-dependent Hamiltonian HM{t) in- 
volves intermittently switched random fields that mimic re- 
peated QND measurements. 

The incoherent S-I cross-polarization transfer that we 
here discuss is then determined by the interplay between 
"free" evolution and measurement effects, as follows: 
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(a) Free -evolution: This will be governed by the time- 
independent terms in Eq. ([T]l- In an interaction picture, 
i.e., in a "doubly" rotating frame with frequencies and 
W/, Hgi will have contributions from both rotating-wave 
{RW or flip-flop) terms S+I^, oscillating as e"("^""^\ 
and co-rotating {CR or flip-flip) terms S^I^ oscillating 
as Q^^i'^s+'^i) ^ and by their respective Hermitian conjugate 
terms. It is clear that the short-time initial evolution will 
be dominated by the rapidly oscillating CR terms, and the 
long-time evolution by their RW counterparts. The energy 
transfer from I to S due to CR and RW terms is then 
governed by the respective population transfer coefficients 

H PcR = 7^Ti£^' = "^^^'^ ^ 

is the effective S — I interaction. At resonance {ijJs = ^i)^ 
Prw = 1> causing a complete exchange of polarization at 
t ~ nvr/ J (n = 1, 2, ...). In such situations, and given that 
usually LOs, ojj > J, one can ignore the fast-oscillating 
CR terms and obtain the dynamics using the RW terms 
only lfl3ll . By contrast, under strongly mismatched condi- 
tions, i.e., ujs ^ oji or vice-versa, Prw ~ Pcb.^ the 
dynamics is equally dominated by the CR and RW terms. 
For such large detuning, the i?5j— driven transfer of polar- 
ization between the S and / spins is inhibited: The polar- 
izations of all spins are then locked at their initial values, 

SLsPcR ~ J/\ujs+uji\,Prw ~ J/\^s-^i\ ^ 1- While 
the presented results are focused on the given Hsi, they are 
general for Hamiltonians that contain RW and CR terms. 

(b) Projective measurements: These will be imparted by 
brief interactions described by H]^j(t) lfl4ll . Each such 
nonselective (unread) projective measurement (l^ erases 
the off-diagonal (correlations) terms that may have arisen 
in the joint S + 1 density matrix. This is equivalent to sub- 
jecting the system to a brief strong dephasing. Although 
the respective eigenstates of the system and the bath re- 
main unchanged during these measurements, their corre- 
lation energy (Hpr) changes drastically, affecting subse- 
quent evolution iflOh . Here, we mimic such projections 
onto the system's energy eigenbasis by an NMR "quantum 
simulator", i.e., spatially-random magnetic field-gradients 
that changes over time (see below). 

The polarization exchange between the S, I spins is 
dramatically altered in the presence of repeated projective 
measurements at times nr, (n = 0, 1, 2, • • • )• To ap- 
preciate this we consider initially uncorrelated equilibrium 
states ps ® pn O ■■■Pin-, products of the 2x2 den- 
sity matrices of the S*— system and each of the N /—bath 
spins that are diagonal in the energy eigenbasis, with pop- 
ulations of the excited spin state being < es(i) ^ 1/2 
(the corresponding polarizations Ps(i) = 1 — "^^8(1))- We 
then find that es(t) oscillates as the weighted sum (over 
all possible /-spin quantum numbers) of S — I oscilla- 
tory exchange probabilities. This function depends on N, 
the bath size, and the anisotropy of the spin ensemble lll4ll 
but primarily on the time between consecutive dephasings: 
(i) At very short times uJs(i)t <^ 1, the S evolution is 



dominated by the fast-oscillating CR terms, so that the 
freely-evolving polarization of the S spin is driven away 
from (1 — 265(0)), causing depolarization of S ("heat- 
ing"): es{t) < es(0). This CR "heating" is amplified 
by the repeated QZE, since CR evolution dominates under 
the QZE condition + W/) r <^ 1 lll4ll : measurements 



or dephasing at intervals Th < ^/V + i^s + W/)^. 
This condition means that highly frequent measurements 
broaden the qubit levels to the extent that they become 
unresolved, equalizing upward and downward transition 
rates regardless of temperature (Fig. [T^: red circles and 
Fig. [ij): lower inset), (ii) By contrast, at longer inter- 
vals, the RW terms increase the polarization (cause "cool- 
ing") of S: es{t) > es(O). Such "cooling", whose con- 
dition is \uJs i uji\t > 1, is amplified by the repeated 
AZE ifioll : measurements or dephasing at intervals ~ 
1 / + {iOs — i-^iY- This comes about since at such in- 
tervals the qubit levels are resolved and downward transi- 
tions (relaxation) dominate at finite temperature (Fig. [T^: 
blue upper triangles and Fig. [TJ): upper inset), (iii) After 
a few measurements (see below) the polarization transfer 
reaches close to the resonant maximum [e/(0)] irrespec- 
tive of the S — I detuning. These time-scales determine a 
resonant-like characteristic that can be exploited as shown 
in Fig. [T] The qubit polarization is then described, within 
the RW domain by Eq. (III.3) in EPAPS IH. 
fcj Quasi-equilibrium: After a few measurements at suit- 
able r's, the polarization approaches its asymptotic value 
and hence the system reaches a quasi equilibrium state (see 
below), with polarization using RW terms only 



,96 



65(0) + 



£/(0) - gs(0) 
2(1-6,(0)) ■ 



(3) 



Depending on the sign of 6/(0) — es(0)' £5^ can be ei- 
ther larger or smaller than £^(0), corresponding to S'-spin 
"cooling" or "heating", as compared to its initial equilib- 
rium value. The value of (1 — 2e'^g) is the largest obtainable 
polarization transfer from the / spins to the S spin, for any 
size of the bath. The transfer achieved by the incoherent 
resonance is thus always greater than 50% of the coherent 
maximum, e/(0)/es(0), and bound by the full coherent 
maximum obtainable under a resonant transfer. 
(d) Reheating: Once e^g is reached, the state of the total 
{S + I) system commutes with the interaction Hamilto- 
nian in the RW approximation, [p, Hsi] ~ 0. This means 
that if no further measurements are performed, the evolu- 
tion of all the spins is almost frozen (Fig. [1^: green lower 
triangles). Yet, in a finite bath as measurements continue to 
be performed, the deviations from Eq. (O due to the CR 
terms, gradually "re-heat" (depolarize) both the 5 and the 
/spins (Fig. lb: blue upper triangles). Hence, different de- 
sired quasi-equilibrium values of the S'-polarization can be 
obtained depending on N, the bath size, and on the number 
of measurements performed beyond the number needed to 
attain eT. 
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Figure 1: (Color online) Evolution of tlie S-spin polarization with time, (a) The main panel compares the evolution of the 5 spin 
polarization, when interrupted by repeated measurements at intervals of t^ ^^ = 1 ms (blue upper triangles) and the evolution interrupted 
at Th = 0.2ms (red circles), respectively. Also illustrated is the quasi-equilibrium state achieved for the Tc^^ measurements, stopped 
after 8 ms and followed by free evolution at later times (green lower triangles). The free evolution of the S spin is shown with black 
squares. The inset zooms the dynamics for short times. These experimental plots are compared with the theoretical curves (dashed lines) 
obtained by exact diagonalization of the Hamiltonian in Eq. (Ujl for the experimental parameters above, (b) Schematic representation 
of the QZE and AZE in thermalized qubits. The white line is the quantum dynamics of the 5* spin steered by 71 = 20 measurements 
for varying the time interval r. It evidences the predicted amplifications compared with the free evolution (black line). For short times 
(QZE regime), the levels are unresolved so that their transition rates are equal (lower inset), while for long times (AZE regime) they 
are resolved so that the downward transitions dominate (upper inset). 



Results. — The foregoing theoretical predictions which 
hold for any size of the bath, were tested by hquid 
state NMR simulators of QND measurement on ^'^C- 
methyliodide (CH3I) dissolved in CDCI3. In liquid-state 
NMR experiments on the methyl group (CH3) a ^'^C spin 
(5) is J-coupled to a finite "bath" of = 3 equiv- 
alent spins (/) which interact with the S spin but 
not with each other. The quasi-equilibrium value of po- 
larization obtained for = 3 is e''g = 65(0) + 
i {(e,(0) - es{0))[l + e,(0)(l - e,(0))]}. The Hamil- 
tonian in Eq. ([T]) was reproduced by applying two radio 
frequency (RF) fields perpendicular to the static field, on 
resonance with the respective / and S spins. In the dou- 
bly rotating frame of the RF fields, precessing with the re- 
spective Zeeman frequencies of the spins we then obtain 
Hamiltonian ([B where the z axis is given by the RF fields 
direction and the frequencies Us and uji determined by the 
strength of the respective RF fields (see lfl4ll for details). 

It is essential that both Ti and T2 relaxation times are 
much longer than the time scales used for these quantum 
simulations. 

To mimic the effects of projective measurements, we re- 
lied on the use of pulsed magnetic field gradients. Field 
gradients effectively increase the decoherence rate for cor- 
relations in a plane perpendicular to the gradient's direc- 
tion, which is along the main Bq field axis ifisll . Yet, it can 
be shown that repeated application of gradients gives rise 
to unwarranted temporal correlations in the coherent evo- 
lution of the system Iil4il . In order to realize a sequence of 
dephasing operations, we resorted to an alternative method 
that truly simulates projective measurements by applying 
at regular intervals new, random, values of the field gradi- 



ents lfl4ll . In an ensemble average, the correlations are then 
effectively erased. 

The initial conditions are [l/2-e/(0)]/[l/2-es(0)] ~ 
4 with the excitations [1/2 — e/(s)(0)] <C 1. By choos- 
ing the time intervals between measurements, to corre- 
spond to the QZE or AZE regimes, depolarization ("heat- 
ing") or, respectively, polarization ("cooling") effects pre- 
dicted by theory were indeed observed under non-matching 
uji ~ 250Hz and cos ~ 420Hz (Fig. [TJi) where the 
choice \ujs — ujj\, Us, uJi ^ J = IbQHz enhances the 
predicted effects. Additional measurements performed af- 
ter attaining the maximum polarization caused re -heating 
of the S spin by CR terms, as theoretically predicted (blue 
upper triangles). Finally, by stopping the measurements 
after maximizing the polarization transfer, we clearly ob- 
served the expected quasi-equilibrium behavior of the S- 
spin (^^C) polarization (green lower triangles). Its value 
agrees well with the theoretically estimated (1 — 2eg) ~ 
2.9 [1 — 265(0)] considering only the RW term. Its slow 
decay is a consequence of non-ideal pulses in the im- 
plementation of the projective measurements. Excellent 
agreement is evident between experimental results and nu- 
merical simulations without any fitting parameters. Fig. [TJ^ 
shows the experimental S-polarization steered by n = 20 
measurements as a function of their time interval r. Heat- 
ing (purity loss) or cooling (purity increase) are seen to 
depend on r as predicted. 

To further explore these incoherent polarization trans- 
fer effects, larger de tunings \ijJs — ^ J i^s ~ 
3.5kHz and Ui ~ 2.6kHz) and fields ujs,ujj ^ J were 
probed and compared with the Hartmann-Hahn resonant 
1 1311 NMR transfer (Fig. |2]l. In these experiments, only 
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Figure 2: (Color online) Time evolution of the system polarization under matching (ujs — ijJj) and off-matching {u)s / loj) conditions, 
(a). The free evolution (black solid line) of the system (S) spin-polarization, and its evolution interrupted by measurements at time 
intervals r of 346/iS (orange circles) and 692/iS (green upper-triangles) are shown for off-resonant fields with high detuning. The given 
r values were chossen to optimize the transfer. The quasi-equilibrium state corresponding to measurements (at time intervals of 692/is) 
stopped after 31.14 ms (n = 45), followed by the free evolution for the later times is marked by blue lower triangles, (b). For resonant 
RF fields, we have plotted the free evolution (black line) of S spin-polarization, evolution interrupted by measurements at time intervals 
1.82ms (green upper-triangle). The quasi-equilibrium state corresponding to measurements (at time intervals of 1.82ms) stopped after 
14.56 ms {n — 8), followed by free evolution at later times, is marked by blue lower triangles. The maximal polarization transfer 
attained by resonant fields (black line) is almost the same as that achieved and maintained (blue lower triangles) by measurements for 
all later times. 



the / spins were initially excited while the S polarization 
was completely erased [es(0) = 1/2] lfl4ll . so that the ac- 
tual polarization transferred from I to S spin could be de- 
termined. Under these off-matching conditions, Prw and 
Pqr were much lower than their resonant values, and the 
transfer arising from the "free" (uninterrupted) dynamics 
was negligible (Fig. |2^-solid line). By contrast, the / — t- 5" 
polarization transfer achieved by repeated projection mea- 
surements was unequivocally evidenced (Fig. ^a)) to be 
close to the maximum achievable under resonant condi- 
tions (Fig.Etb)). 

Discussion. — The theoretical and experimental data 
shown in Figs. [T] and |2] clearly demonstrate for the first 
time the connection between the fundamental Zeno and 
anti-Zeno effects in frequently measured or dephased open 
systems and their purity loss or increase, respectively. 
Measurement (dephasing)-induced transfer of polarization 
has been shown to be almost as effective as coherent on- 
resonance transfer, even if the measured system and bath 
are under presumably unfavorable, off-resonant conditions. 
We have further demostrated the ability to steer the system 
into a quasi-equilibrium, which is maintained when further 
measurements are stopped. 

In terms of their practical use, we envisage potential ap- 
plications of this non-unitary polarization transfer protocol 
for qubit purification, required at the initialization stage of 
quantum information processing ||l^. The possibility of 
increasing the polarization transfer from the pure / spins 
to the impure S spins even under unavoidable off-resonant 
conditions could be useful for algorithmic cooling ifivll . 
Hence, the applied and fundamental aspects of the incoher- 



ent resonances introduced here merit further exploration. 
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Supplementary information: Zeno and anti-Zeno polarization control of spin-ensembles by 

induced dephasing 



I. QUBIT DYNAMICS WITHIN THE ROTATING WAVE APPROXIMATION 



Let us consider the Hamiltonian of Eq. (1) of the main article within the rotating wave approximation, i.e. 



(1.1) 



The total Hamiltonian has a 2 x 2 block-diagonal structure, each block having a definite value of S^+P. Its time-evolution 
operator has the form 



N/2 I r 

7=0 Mi = -I I 



(7 + T^CT-^ 



Mi 



Vi, 



(1.2) 



where Mj is the magnetic quantum number of /, Ja/^ = J^y {I — Mj){I + Mj + 1), A = (ujs — W/)/2 and 



Jfj^ + A^. (t's are the Pauli operators in the basis of the 2x2 block diagonals of Him,-- The projection operator Vi 
corresponds to each bath-spin sector. The phase (pMi corresponding to each block does not contribute to the dynamics, if 
the initial state is diagonal in the 5^ + P basis. 

The initial density matrix in the basis of | i; Mj), | — |, Mj + 1) is given by 



and the diagonal element evolutions are 



ee 



I. 99 , 



1 



P 

'-'Mi 



sin^ Mi t 



e5(0)(l-e,(0)) 

e,(0)(l-es(0)) 



ee 1 , \gg 
k99 i Jhi ■ 2n ^1 

PMi +^Sm ilMitpMi 



(1.3) 



(1.4) 



Mi 



Complete exchange of polarization between / to 5", i.e. es(r) = e/(0),is only possible at times r . This coherent exchange 
is destroyed in presence of off-resonant fields (A ^ 0). In general, the coherent exchange is controlled by the transfer 
coefficients 



P+_(n) 



and the spread of the incommensurate frequencies O Mi ■ 



, l<n=J{I-Mi){I + Mi + l)<VN, 



(1.5) 



A. Counter-rotating terms 



When the interaction Hamiltonian has counter-rotating (CR) terms only, i.e., 

Hen = oosS' + E + ? E [S^^>^ + S'^^] ' 



(1.6) 



then + is no longer a conserved quantity. Yet, the Hamiltonian can be expressed in a block-diagonal form in the 
basis of I i; M/), | — |, Mj — 1) . Similarly to the above analysis, the time-evolution of the total system exchange energy 
between S and I determined by the factor 



(1.7) 
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II. REPEATED MEASUREMENTS 



A measurement on the system which erases its off-diagonal elements while keeping populations unchanged is a nearly 
ideal quantum-non-demolition (QND) measurement. This measurement projects the system's state onto its energy eigen- 
basis. In the present case this is also equivalent to projecting also the /-spins on their (total) basis. For the RW and 
CR interaction dynamics, the projective measurement can be performed individually in each subspace of + I^' . By 

Xm, 

1 - XMi 

elements are erased by each projective measurement, if n of them are performed at time-intervals r, then 



rewriting the initial density matrix in Eq. (11.31) as = Tr[pj\/j] 



and accounting that off-diagonal 



.M.(i = nr) = /rM,(r)xM. + /2M.(r) J] /rM,(r), 



m=0 



2.n 



- sin^ ^Mjt and f2(t) = sin^ ^Mit. In the limit of n — )• oo, xi\ij (t) 



where fi{t) = 1 



the density matrix p'/^ = pMj{t — oo) = Tr[pA/J 



(II. 1) 



h 



Hence the total density matrix (qubit -i- bath) commutes 



with the Hamiltonian A similar equilibrium state can be found for the Hqr- 



ni. STEADY-STATE VALUES OF THE QUBIT POPULATIONS 

The time-dependent population of the S spin can be found by tracing over the bath degrees of freedom obtaining 

where F{t) = j,^^ ^lij sin^ ^M,t and < F{t) < l.The weight function W/j^ is given by 



(III.l) 



(ni.2) 



where A/ = (jv/2+j)!('jv/2-/)! n/ 2+1+1 • performing repeated measurements as in the above section, the final S'-spin 
polarization within RW approximation is given by 



1 - es(0) 



(III.3) 



where 



Mr) 



2{eji0) - esm fj^^ 
26,(0)(l-e,(0))(l-r(T))]. 



E WL, [ {6.(0)(1 - 6.(0)) + 65(0)(1 - 6,(0))} (1 - fir)) J2 r{r) 



m=0 



(111.4) 



In the above equation f{T)= cos 2Q]\ij r. After few measurements /" (r ) becomes negligibly small and J2m=i3 /™ (''") 
jzj^- Then, a steady state value is attained for any bath size N, given by 



esiO) + 



6,(0) - 65(0) 
2(1-6/(0)) 



1_ J^A,6,(0)^+^(l-6,(0))^-^ 



The equilibrium value for, = 3, can be evaluated from the above equation, giving 



1 



65(0) + - {(6,(0) - esiOm + 6,(0)(1 - 6,(0))]} 



(III.5) 



(ni.6) 



Since < 65(0), 6,(0) < 1, there will always be a gain in polarization 61*^ > 65(0), if 6,(0) > 65(0). If the interaction 
Hamiltonian has counter-rotating terms only, the projective equilibrium value is given by 



1 



[ef]cR = esiO) + - {(1 - 6,(0) - 65(0))[1 + 6,(0)(1 - 6,(0))]} . 



(111.7) 
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Thus the counter-rotating terms take es(0) close to 1 — e/(0). 



Large N limit 



For large one can attain the maximum achievable polarization transfer from the bath spins to 
the qubit (using RW terms only). This saturation value, e^^ ~ ^^•(O) + ^^r^f^, is reached for any 
value of the magnetic fields, us and uj. 
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Figure 3: NMR pulse sequence, (a) Both kinds of spins are excited for the initial condition (Fig. 1, main text), (b) Based on the phase 
cycling of pulses and the acquisition process, only the I spins are effectively excited as initial conditions of this variant (Fig. 2, main 
text). The time length of the 90° pulses were 8 for the / spins (^H) and 15 for S (^'^C). We used for the random intensities of the 
gradient values in the range \ABz | < 30 G cm, where the length of the sample is approximately 1 cm. This leads to a variation of the 
precession frequency around the static field of approximately 125 kHz for and 32 kHz for ^'^C. The time duration for the gradients 
is Tm — 100 ^s, satisfying the condition Aoj/jTm — ■js^Bzhrm 3> 2it. 



IV. EXPERIMENTAL DESCRIPTION OF THE NMR SIMULATOR 



A. Hamiltonian realization 



The experiment described in the main text involves one spin S and three bath spins / in the presence of a static magnetic 
field Bq in the z direction. Considering that the differences of the Larmor frequencies are far greater than the J-coupling 
between the spins, the effective static Hamiltonian is given by 

H = ojioP + cosoS' + JPS% (IV. 1) 

where P = I^, the sum being over bath spins. When radio-frequency (RF) magnetic fields Bu and Bis oscillating 
at frequencies ujjq and ojso are applied, the Hamiltonian in a doubly rotating frame, processing with the respective Larmor 
frequencies, becomes O 



H = ujiP +ojsSy + JPS\ 



(IV.2) 
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where = jiBu and ^5 = 75S15. 

We assume that the direction of the RF fields is — y in this frame. By making a change of variables for the spatial 
directions — )■ z, z — )• x) this Hamiltonian becomes equivalent to the one given in Eqs. (1) and (2) of the main text, 

H = uiP + uJsS' + jrS\ (IV.3) 

We would like to note that the presence of both RW and CR terms in the above Hamiltonian is a direct consequence of 
the anisotropic I^S'' coupling. The I ■ S ox S^I^ + S^P coupling, as found in many natural spin systems, would result 
only in AZE cooling and quasi-steady state behavior discussed in the main text. The QZE heating and the reheating effects 
would be absent, because of the cancellation of the CR terms. 



B. Pulse sequence with random gradients 

Figure [3] is a schematic representation of the NMR pulse sequence used for simulating the measurement- induced quan- 
tum dynamics described in the main text. The preparation of the initial state, shown in Fig. |3] consists of the following 
steps: 

(i) . 90° RF pulses applied along the x axis jsl], exactly on-resonance with the respective Larmor frequencies of both 
spin species. These pulses rotate the magnetizations from the initial z direction to the —y direction. This is followed by 
the application of on-resonance RF fields pointing to the —y direction on both kinds of spins, and the Hamiltonian iYV.3i 
is obtained during time r. 

(ii) . In order to simulate the projective measurements, we rotate the polarization back to the z axis by 90° pulses on the 
—X axis. At this point a static field gradient is applied during time t„i. To simulate a projective measurement, the gradient 
intensity is chosen randomly (see section IVll. 

Thus the Hamiltonian generated by the field gradients in a slide of the sample at a position z on the axis parallel to the 
gradient direction is given by Hlj{t) oc J2n ils^B'zZSz + jiAB^zI^) 9 (t — nr) 9 (nr + t,„ — t) , where Ai3" is 
the gradient strength at the step n. After each r„j time, the ensemble average of the different portion of the sample preserves 
the populations in the density matrix but erases the quantum correlations between the / and S spins. Additionally, the 
non-diagonal elements of the density matrix within the blocks of constant total spin of /, the so-called zero quantum 
coherences, are not erased, due to the full chemical equivalence of the / spins. 

(iii) . Following this projection the polarizations are returned to the —y direction by 90° pulses on the x axis, and the 
entire cycle is repeated as schematized in Fig. |3] After n cycles, we realize the acquisition of the Sy magnetization as if n 
projective measurements were performed on S. As the coherences between / spins are not modified by the gradients, the 
simulated measurements mimic measurements only on the S spin. 

(iv) . Alternatively one can prepare an initial excitation on the / spins only. For this purpose, during the preparation 
time, only the / spins are rotated to the —y direction, keeping the S magnetization on the z axis. By phase cycling we 
cancel out the initial contribution of the S spins (i.e., €^(0) = 1 /2). 

C. Pulse sequence without random gradients 

In order to compare these evolutions with the free-evolution dynamics, we repeated the experiment with the same 
pulse sequences but with null gradients. This procedure gives the same time scales for the relaxation and for other non- 
ideal features of the pulses, but does not modify the free dynamics between measurements. Decoherence is manifest 
by the reduction of the oscillation amplitude as time goes on. However, the first oscillations are comparable for both 
procedures. This provides the evidence that there is no polarization transfer during the tt/2 pulses. The polarization 
control by simulated measurements discussed in the paper is only achieved when the random gradients are introduced to 
dephase of the coherences (correlations) between S and / spins (Eq. |V.2| |. 

V. GRADIENT-DRIVEN PROJECTIVE MEASUREMENTS IN NMR 

The experimental simulations of repeated projective measurements were achieved with random gradient fields. A brief 
analysis of this equivalence is given below. 

We consider a S spin in presence of a magnetic field on the y axis whose Zeeman frequency is Its initial state is 
Ps = exp{—/3u}QsS~)/Z ~ (1 — I3ujqsS') jZ in the high temperature limit, a typical condition in NMR experiments 
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Its evolution is given by 

p (r) = 5^ cos (wst) + 5^ sin (cjst) , (V.l) 

where we do not write the constant evolution of the identity operator and the factor —f3ujQs/2. By performing projective 
measurements at time intervals r on the energy eigenbasis of 5^ , S"^ sin (w^r) is erased. After n measurements p„ (nr ) = 
cos" (wgr) . Applying a field gradient on the z axis at time intervals r during a time t,„ to simulate the measurements, 
a slice of the sample at position z is affected by the field AB^z, where Ai?^ is the gradient value. The resulting evolution 
of the slice's density matrix (IV. lb is (r) = cos (w^r) + [5^ cos [Aw^Tm] + sin (Aw^Tm)] sin (w^r) , where 

Aojz = ^s^B^z. Its ensemble average within the sample length h is then p (r) = /_''(/2 (''") ^ ~ (^st) ; 

assuming that Aw/jTm = ^s^B^hr^ S> 27r. 

While this result represents a projective process for n = 1, the next step n = 2 is no longer the outcome of two 
consecutive projections. By calculating the density matrix of the portion z of the sample at this step and considering the 
ensemble average 

rh/2 nh/2 

p (2t) = / pz (2t) — = S^' cos^ {(jJst)+S^ sin {(jJst) / [cos (ojst) cos^ (Aw^Tm) — sin^ (Ac^j^Tm)] — . 

(V.2) 

While the first term on the rhs is equivalent to the one obtained by two projective measurements, but the second term is 
non-vanishing. In general, correlations of multiple steps, of the form cos™ (Aaj^Tm) sin' (Aa;.ri„) with / and m even 
numbers, survive the ensemble average, leading to non vanishing deviations from projective measurement. Therefore, a 
field gradient by itself does not simulate repeated projective measurements. By changing randomly the gradient intensities 
at consecutive steps, we destroy correlations among different steps upon taking the ensemble average. Thus only the term 
S"^ cos" (wst) survive the ensemble average, properly mimicking the projective measurements. 

Similar results are obtained once considered the full S-I system. The Hamiltonian Hli{t) described in the preceding 
section dephases Sx.y and I^.y terms in the density matrix when the ensemble average on the position z is taken. An 
ideal projective measurement is performed when all spins along z experiences a random field. A single gradient gives a 
particular z-dependent phase to spins along the sample. Only when this gradient is randomized further many times such 
that the sum of the phases accumulated under each gradient vanishes for each spin, one can realize an ideal projective 
measurement as showed above. However, because the non-diagonal elements of the density matrix within the blocks of 
constant total spin of / are not erased, due to the full chemical equivalence of the / spins, the simulated measurements 
mimic measurements only on the S spin. 
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